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The quality of supernova data will dramatically increase in the next few years by new experiments 
that will add high-redshift supernova to the currently known ones. In order to use this new data to 
discriminate between different dark energy models, the statefinder diagnostic was suggested pj and 
investigated by Alam et ai.0| in the light of the proposed SuperNova Acceleration Probe (SNAP) 
satellite. By making use of the same procedure presented by these authors, we compare their 
analyzes with ours, which shows a more realistic supernovae redshift distribution and do not assume 
that the intercept is known. We also analyzed the behavior of the statefinder pair {r,s} and the 
alternative pair {s,q} in the presence of offset errors. 
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Introduction 

Recent observations from type la supernovae measure- 
ments, cosmic microwave background radiation and grav- 
itational clustering suggest the expansion of the universe 
is accelerated. 

In order to explain this cosmic acceleration a form 
of negative-pressure matter called dark energy was sug- 
gested. The simplest and most popular candidate is 
Einstein's cosmological constant. Many others candi- 
dates for dark energy have been proposed, including 
scalar fields with a time dependent equation of state, 
quartessence, modified gravity, branes, etc. Confronta- 
tion between these models and currently observational 
data doesn't say much mainly because most of them 
have ACDM as a limiting case in the redshift range al- 
ready observed. The SNAP (SuperNovae Acceleration 
Probe) satellite is expected to observe ~ 2000 supernovae 
per year with redshift up to z — 1.7. To differentiate 
models using the new available data, Sahni et al. 0] in- 
troduced the statefinder diagnostic, that is based on the 
dimensionless parameters {r,s}, which are constructed 
with the scale factor and its time derivatives. 

In this work we applied the statefinder to a SNAP-like 
supernovae distribution and analyzed its behavior in the 
presence of systematic and random systematic, beyond 
statistical errors. 



energy with an equation of state in the form p x — p x {p x ). 
We also take c = 1 and consider a flat universe (k = 0). 



The focus of our discussion will be in the four models 
listed below: 

1. Cosmological Constant (w x — p x /p x — — 1) 

The cosmological constant model represents a con- 
stant energy density. In this model, the Hubble 
parameter has the form: 



H(z) = H [ Ono(l + zf + 1 - n m0 p 



(3) 



2. Quiessence ( —1/3 > w x = p x /p x = cte > —1) 
This is the next simplest example of dark energy 
model, and gives rise to a Hubble parameter like: 

H(z) = H [ Q m0 (l + zf + njro(l + z) 3 ^]i (4) 

3. Quintessence ( w = w ( t ) ) 

Representing a self-interacting scalar field mini- 
mally coupled to gravity. In this model, we have: 



P4> = \j>{zf + V{ct>{z)) 
P<t> = -V{4>(z)) 



(5) 
(6) 



I. DARK ENERGY MODELS 

Assuming a Friedman-Robertson- Walker (FRW) met 
ric, the Einstein's equations reduce to: 



H 2 = 
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where a is the scale factor of the FRW metric, H is the 
Hubble parameter, and the sum is over all the compo- 
nents present in the scenario in study. 

In the following we assume that the matter content of 
the universe is given by dark matter (p m = 0) and dark 



We shall focus on a special kind of quintessence 
model that has a tracker like solution, with the 
following potential: V((p) — (f)(z)~ a (a > 0). In 
this case, it can be shown that the present energy 
density of the dark energy is almost independent of 
initial conditions. 

4. Chaplygin Gas 

A different kind of solution is provided by the Chap- 
lygin gas model. In this model, the Hubble param- 
eter takes the form: 



H(z) = H a 



zf + 



a, 



(7) 
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Table I: Redshift distribution. The value of z represent the upper edge of each bin |J 



z 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


N(z) 





35 


64 


95 


124 


150 


171 


183 


179 


170 


155 


142 


130 


119 


107 


94 


80 



where 
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The Chaplygin gas behaves like a cosmological con- 
stant for small z (late times) and like pressureless 
dust for large z (early times). 

It is important to note that for all models presented 
previously, the luminosity distance is : 



Dl{z) 
1 + z 



dz' 



(9) 



with H(z) given by the model dependent expressions pre- 
sented before. 



II. THE STATEFINDER 

The properties of dark energy, as we have seen, are 
very model dependent. In order to differentiate be- 
tween the presented models, Sanhi et al.[l], proposed 
the statefmder diagnostic. The parameters, {r,s}, are 
a complement to the already known deceleration param- 
eter, and help the discrimination when the later contains 
degeneracies. By definition: 
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III. DARK ENERGY FROM SNAP DATA 

Type la supernovae are considered standard candles, 
used to map the expansion history, and its observations 
lead to the behavior of the scale factor with time. In 
order to study the data in a model independent way, 
we use a parametrization for the dark energy density, 
presented by Sahni et al. We express the energy 

density as a power series up to second order in z: pde — 
Pco{Ax + A 2 x + A 3 x 2 ), where x — 1 + z. For this Ansatz, 
the Hubble parameter takes the form: 



equation ljl"3j) together with equation JHJ provides an ex- 
pression for the luminosity distance, which we shall in- 
vestigate, using the statefmder parameters, in the light 
of a SNAP-like experiment simulation. 

To simulate the data, we use a binned approach for a 
SNAP distribution shown in the Table I. We also include 
300 supernovae in the first bin. These low redshift super- 
novae are expected from the SNFactory (Nearby Super- 
novae Factory) experiment and are important in reducing 
the systematic errors. The SNFactory proposal is to pro- 
vide data to calibrate high redshift experiments, like the 
SNAP, and then reduce the errors involved. 

The luminosity distance of equation © is measured in 
terms of the apparent magnitude, which can be written 
as: 

m(z) — 5 \ogd(z) + 

[M + 25 - 5 log (H / (100km/ s/Mpc))]{U) 

where M is the absolute magnitude of the supernova and 
the expression in brackets is called intercept. 

Following what was done by Kim et al. jj| , we consider 
a random irreducible systematic error of 0.04* (z me d/1.7) 
(here z me d is the redshift in the middle of each bin), 
added in quadrature to a constant statistical error of 
0.15mag and study the behavior of the statefmder pa- 
rameter with this synthetic data. We performed a Monte 
Carlo simulation considering the intercept exactly known 
and totally unknown. As a second step, we study the 
situation where offset errors are present, and its conse- 
quences in the statefmder. 



IV. RESULTS 

In the figures to El we present the results from simu- 
lated data. According to SNAP's specifications, we gen- 
erated 500 data sets having ACDM as a fiducial model. 
For each of this experiments we calculated the best fit- 
ting parameters A\ and A 2 for equation 11 -tt and recon- 
structed the statefinders r(z) e s(z). The figures below 
show the mean value of the parameters, which were cal- 
culated as: 



H(x) = H Q (n m0 x 3 + Ai + A 2 x + A 3 x 2 )i 



(13) 
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Figure 1: Shows < r > as a function of redshift. The blue contour 
represent la and the red contour 2a confidence levels, in the pres- 
ence of statistical errors with a known intercept. The line < r >= 1 
represents the ACDM fiducial model. The dashed lines above the 
ACDM are Chaplygin gas with parameters k = 1 and k = 2. The 
dotted lines below it are quiessence models with w = —0.6 and 
w = —0.8, and the dotted-dashed lines are quintessence models 
with a = 2 and a = 4. 
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Figure 2: Analog to figure Q but here in the presence of random 
systematic and statistical errors, with an unknown intercept. 
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Figure 4: Shows < s > as a function of redshift in the presence of 
statistical error. The colored contours, dotted, dot-dashed and full 
lines represent the same models as in figure Q 



l 

0.75 
0.5 
0.25 


-0.25 
-0. 5 
-0. 75 



o= 4 . 




QUINTESSENCE 

— ■ - — *i= -tr: 6 ; 
















CHAPLYGIN GAS 

k= 2 



0.25 0.5 0.75 1 1.25 1.5 



Figure 5: Analog to figure Bvl but here in the presence of random 
systematic and statistical errors, with an unknown intercept. 
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Figure 6: Shows < s > as a function of redshift. Again, the 
blue and red contours are fcrand 2a confidence levels, in this we 
considered a known intercept in the presence of statistical error 
and a systematical errors of +0.03macj.The dotted, dot-dashed and 
dashed lines are the same as in figure HI 



Figure 3: Shows < r > as a function of redshift. Again, the 
blue and red contours are la and 2a confidence levels, in these 
we considered a known intercept in the presence of statistical error 
and a systematic error of +0.03mag.The dotted, dot-dashed and 
dashed lines are the same as in figure 



Next, we investigated the integrated averaged of the 
cosmological parameters, as suggested by Alam et al. 
Q. For the cosmological constant model, the parame- 
ters are constant, but for many dark energy models the 
statefinder evolves (as it is clear in the figures presented 
previously), the integration of this quantities may then 
reduce the noise in the original data. The integration 
was done as follows: 
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where z max = 1.7 . The expressions for r, s and q were 
calculated using equations El Ell and the first part of 
equation El for each experiment. The 500 points ob- 
tained are plotted in figures HIH to 1151 
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Figure 7: Shows < q > as a function of redshift. The blue contour 
represent la and the red contour 2a confidence levels. The dashed 
lines are Chaplygin gas with parameters k = 1 and k = 2. The 
dotted-dashed lines are quintessence models with a = 2 and « = 4. 




-0. 75 



0.25 0.5 0.75 1 1.25 1.5 



Figure 8: Analog to figure but here in the presence of random 
systematic and statistical errors, with an unknown intercept. 




-0. 75 



0.25 0.5 0.75 1 1.25 1.5 



Figure 9: Shows < q > as a function of redshift. Again, the 
blue and red contours are la and 2a confidence levels, in these we 
considered a known intercept in the presence of statistical error and 
a systematical error of +0.03mag.The dot-dashed and dashed lines 
are the same as in figure 171 



V. DISCUSSION 



Figure 10: Shows the variation of f with s , in the presence of sta- 
tistical error when the intercept is known. The red dot is the ACDM 
fiducial model. The blue dots above the ACDM are quiessence mod- 
els with w = 0.0, —0.3, (—0.5), —0.7 and 0.9 respectively, from top 
to bottom. The green dots below ACDM are Chapligyn gas with 
k = (3/7), 1, 2, 6, 15, from bottom to top. 
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Figure 11: This is analogous to figure Hul although here we consid- 
ered an intercept not known and included statistical and random 
systematic errors. 
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Figure 12: Show the variation of f with s when the intercept is 
known, in the presence of statistical error and systematic error of 
+0.03mag. The red, blue and green dots are the same as in figure 



Our results suggest that the statefinder is a good diag- 
nostic for dark energy models, although some care must 
be taken in applying it to data. As expected, the pres- 
ence of random systematic error and an unknown inter- 
cept just added more possible models than those allowed 
by statistical errors only. There is little problem in this, 
once the fiducial model is at least within the 2er contours 
(Figures 00 and El- 

The presence of offset errors had different outcomes. 
For a positive offset, the parameters r and s suffered 



a small reduction in relation to the fiducial model, but 
in different redshift ranges (r for small and s for high 
redshift). A consequence of this arises when we compare 
the integrated averages of the pair {r,s}. In Figure El 
the points are shifted to the negative direction of both 
axes of the ellipse, resulting a data set where the fiducial 
model (red point) is on the edge of the distribution. The 
same kind of shift is observed for a negative offset, but 
in this case to the positive direction of the axes. So, in 
order to use the statefinder as a diagnostic, we must be 
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Figure 13: Shows the variation of s with q, in the presence of statis- 
tical error when the intercept is known. The red dot is the ACDM 
fiducial model. The blue dots above the ACDM are quiessence 
models with w = 0.0, —0.3, (—0.5), —0.7 and 0.9 respectively, from 
top to bottom. The green dots below ACDM are Chapligyn gas 
with k = (3/7), 1, 2, 6, 15, from bottom to top. 



able to control the offset error below 0.03mag. 



It is interesting to observe the behavior of the decel- 
eration parameter when systematic errors were involved. 
It does present a very small reduction (positive offset) at 
high redshift, but it is irrelevant in front of that suffered 
by r or s. Comparing figures an d E3I we could say, as 
suggested by Alam et at 0], that the pair {q,s} is even 
a better diagnostic than {r,s}, once it restricts the area 
of the phase space filled by the data. However, if there is 
a systematic error present, the points will be shifted to 
the negative direction of the q axis only (figure El > let- 
ting the distance between the data and the fiducial model 
bigger than those in figure IT21 
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Figure 14: This is analogous to figure 1131 although here we consid- 
ered an intercept not known and included statistical and random 
systematic errors. 



Therefore, we concluded that the statefinder pair {r,s} 
is a better diagnostic than the pair {s,q}, when the in- 
volved systematic errors are not random. 
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Figure 15: Show the variation of s with q when the intercept 
is known, in the presence of statistical and systematic error of 
+0.03mag. The red, blue and green dots are the same as in figure 
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